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Abstract. The electromagnetic radiative correction to the cross section of the vector meson electroproduc¬ 
tion is calculated. Explicit covariant formulae for the observed cross section are obtained. The dependence 
of the radiative correction on the experimental resolution and on the inelasticity cut is discussed. FOR¬ 
TRAN code DIFFRAD based on both exact (ultrarelativistic) and approximate sets of the formulae for 
the radiative correction to the cross section is presented. The detailed numerical analysis for kinematical 
conditions of the recent experiments on the diffractive electroproduction of vector mesons is given. 


1 Introduction 

The measurement of the cross section of the exclusive vec¬ 
tor meson electroproduction can provide information on 
the hadronic component of the photon and on nature of 
diffraction. During several years the diffractive production 
of the vector meson has been the subject of the muonpro- 
duction @,11 and electroproduction experiments. 

Data analysis of these experiments is considerably affected 
by the QED radiative effects. At practice the radiative 
corrections (RC) to the processes of electroproduction are 
taken into account using codes originally developed for the 
inclusive case (see ref.0, for example). 

The purpose of this paper is to calculate the electro¬ 
magnetic correction to experimentally observed cross sec¬ 
tions for the kinematics of a fixed target and collider ex¬ 
periments directly. The Feynman diagrams necessary to 
calculate RC are presented on Figj^. 

In order to calculate exactly the QED RC to the cross 
section of vector meson production, the method offered in 
ref.|^ is used. By exact formulae is meant the expressions 
for the lowest order RC obtained without any approxima¬ 
tions but ultrarelativistic: the lepton mass m is considered 
to be small. In Section ^ the kinematics of the radiative 
and non-radiative processes and exact formulae for the 
lowest order RC are obtained. In Section || the analytical 
results are visualized by the construction of the approx¬ 
imate formulae for cases interesting in practice. The nu¬ 
merical results are given in Section 0 A brief discussion 
and conclusions are given in Section g. 

2 Exact formulae for the lowest order 
correction 

Seven kinematical variables are necessary to describe the 
radiative process of the diffractive vector meson produc- 



d) e) 


Fig. 1. Feynman diagrams contributing to the Born and the 
next order cross sections. The letters denote the four-momenta 
of corresponding particles. 

tion (Figure 03, c). Four of them are the same as for non- 
radiated case: usual scaling variables x and y, the negative 
square of momentum transfered from the virtual photon 
to the proton t = {q — Ph)^ and the angle (j^h between 
the scattering (ki,k 2 ) and production (q,Ph) planes in 
the laboratory frame. The squared virtual photon mo¬ 
mentum = —{ki — k 2 )^ and the invariant mass of 
the initial proton and the virtual photon = {p + q)^ 
are often used instead of x and y. The kinematics of a 
real photon is described by three additional variables: the 
inelasticity v = (M is the proton mass and 
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= {p' + kY is the squared invariant mass of the system 
of unobserved particles), t = kq/kp and angle 4>k between 
planes (ki,k 2 ) and (q, k). 

We consider RC to three and four dimensional cross 
sections tr = da/dxdydtd(j)h and d = da/dxdydt. They 
are related as 

27r 

a = y* d(t)h (7 . (1) 

0 

The four differential Born cross section can be presented 
in the form 

^ 4^ + 2(1 - y - \y^l^) {OL + ar )^ , (2) 


where ut and ctl are differential cross sections of the pho¬ 
toproduction, 7 ^ = and v is the virtual photon 

energy. 

The differential cross section of the radiative process 
has the following form: 


Wh 




where A = p + q — ph, are matrix elements of the 
processes given in Figure |po,c. In order to extract the 
infrared divergence into a separate term we following Q, 
use the identity 


= ctr - <yiR + <yiR = ap + aiR, (4) 

where Ui? is finite for k 0 and am is the infrared diver¬ 
gent part 


(tir = -(5jj^cro = -(^s + SH)ao. (5) 

TT TT 

The quantities Ss and Sh appear after splitting the inte¬ 
gration region over v by the infinitesimal parameter v: 


0 

Vm 

^H = lJdvJ ^Fir 5{{A - kf - M2), 


m2). 


( 6 ) 


where p is an arbitrary parameter of the mass dimension, 
Vm is a maximal inelasticity and 


Fir 


2 2 
m m 

(2fcfci)2 (2 A:/c2)2 


(52 -(- 2m^ 

{2kki){2kk2) 


(7) 


The wiw to calculate the integrals like has been offered 
in ref. [g (see also review [§ ). In our case we have 


Ss = 2 (P,s + I 08 (i™ - 1) + log 4^ + St. 

= 2(/^-l)log^, 

V 


where Im = \og{Q'^/m^). The quantities S' = 2Ak\ = 
S — Q'^ — Vi and X' = 2Ak2 = X + Q'^ — V 2 are calculated 
using Fi ,2 = 2 ( 01,2 + b cos (j)h), where 

- (55, + 2M^Q^)tg), 

ZAq 

02 = T^iQ^SpSt - (X5, - 2M^Q^)tg), 

ZAq 

b = - StSj, - Mhl - 

Aq 

x(5XQ2_m2q4_to2^ji/2^ (9) 

The invariants are defined as 

5 = 2kip, X = 2k2P = (1 - y)S, Q2 = Sxy, 

Sp^x = 5 ± AT, St = Sx + tq = t + 

Aq = 52-k4M2Q2. (10) 

The infrared terms Pjr, parameters p and v as well as 
the squared logarithms containing the mass singularity 
are completely canceled in the sum of (5|j^ with 6v coming 
from a contribution of the vertex function (Figure 0d): 

by = —2(F7ij-t-log—)(lm —1) —(II) 

p 2 2 D 

For this sum we have 


= binf + bvR, 


( 12 ) 


where 


bvR = 


a 


3, „ 1 , 2 2^' 

2 ™ “2- 2 log 


LE 


I - 


Q2M2 

S'X' 


TT 

"6 


bsnf=^{lm-l)l0g^. (13) 

Here the ultrarelativistic expression for 5^ calculated in 
[ pT was used. The higher order corrections can be par- 
tia ly taken into account using a special procedure of ex¬ 
ponentiation of the multiple soft photon radiation. There 
is an uncertainty, what part of byu has to be exponen¬ 
tiated. Within the considered approach Q (I -|- bmf) is 
replaced by exp Sin f ■ 

For the observed cross section of the vector meson elec¬ 
troproduction we obtain 

o-Qbs = croe^*"^(I -I- bvR + byac) + a-p- (14) 


The correction Syac comes from the effects of vacuum 
polarization by leptons and hadrons (Figure^). The ex¬ 
plicit QED formulae for the hrst one can be found in 
ref. . The hadronic contribution is given by a fit coming 
from the data on e“'"e“ —> hadrons p 2 [ . 

The contribution of the infrared hnite part can be writ¬ 
ten in terms of POLRAD 2.0 notation 


q;2j/ 


Tmax 2 3 


7’"-^ 


Q4 


(15) 


( 8 ) 


0 
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where R = v/f, f = 1 + t - and 2MTjnax,min = ± 

Sj = 1 for j = 1 and 6j = 0 otherwise. 

The quantities Oij depend only on the kinematical in¬ 
variants and the integration variables r and (j)k 


to make a cut on the maximal value of this quantity. If 
this cut is not applied the maximal inelasticity is defined 
by kinematics only. Below we give the formulae for Vm in 
terms of the kinematical invariants 


011 = 

012 = 

013 = -2{2F + FdT^), 

021 = A{SX - M^Q^)F°^, 

022 = -FdS^T + Fi+SpS, + 2F2-Sp 

023 = AFM^ + 2FdM^T^ - FdS^T - F^+Sp, 


where 


AQ'^Vm = -b - 

-{Sx-2Q^ +tqf -AM^Q^ (21) 


(16) 


and in terms of kinematical limits on t 


1 


— ^{imax 0 ^mm); 


( 22 ) 


where C behaves for small t as 


Fd 


F F F 

-, J^i+ =-1-, F2± = F 

Z 1 Z 2 Zi Z2 



(17) 


C = 


+ ml 
21T2 




+ 0{t). 


(23) 


F = 1/{2tt^) and = F 2 + - Q^Fd = FmR^- The 
quantities /r = kph/kp = 2{ak+bk cos{(j)h — 4'k)) and zi _2 = 
kki^ 2 /kp = 2(of 2 ~ cos(f>k) include the dependence on 
angles: 


The maximal inelasticity given by kinematics are plotted 
in Figure 


a? = 77^(Q"5 'p + TiSSx + 2M^Q^)), 

ZAq 

“2 = T^iQ^Sp + t{XS, - 2M^Q^)), 

ZAq 

ak = ■^i(2Q^ + TS^)Sm - - 2rM^)t,), 

ZAq 

F = ^{SXQ^ - 

A, 

x(Q 2 + 1-5-^ _r2M2)i/2^ 

bk = - S^Sxtq - MHI - mlXqf/^ 

x{Q^ (18) 

Here Sm = St — v. 

The dependence on the photoproduction cross sections 
is included in Fp. 



F\ — {Sx — R)crT^ F 2 — g _ 

Fi = SxOt, T 2 = 2x{aT + ctl). 


Fig. 2. Maximal inelasticity within kinematical conditions of 
(19) HERMES. Dashed line gives a possible cut. 


The quantities ctt.l have to be calculated for Born kine¬ 
matics, but is calculated in terms of so called true 
kinematics. It means that they have to be calculated for 
the tilde variables 


q 2 = g2 ^ 

-R(1 + t), 

i = t + R(t — /i) (20) 

instead of usual Q^, and t. 

The important point is the dependence of the results 
on the maximal inelasticity Vm- The inelasticity is calcu¬ 
lated in terms of the measured momenta, so it is possible 


3 The approximation 


The structure of the exact formulae for the cross section of 
the hard bremsstrahlung ap obtained in the last section is 
not simple. Here we analyze them under approximation. 
The third integration variable v is not completely photonic 
one. It is calculated from the measured momenta of the 
final lepton and the vector meson. For non-radiated events 
z) = 0. At practice, however the distribution over v has 
a Gaussian form due to the finite resoluti on. Radiative 
effects also affect the distribution (see sect. O). So a cut 
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on the invariant mass of the unobserved system rejects 
events with hard radiated photon and helps to reduce the 
total RC. Below we consider the approximation based on 
assumption that Vm is relatively small or kinematical cut 
on V is used. Besides, for the diffractive scattering the 
squared momentum t transfered to the hadronic system is 
always small, so we construct the approximate formulae 
for RC to the cross section da/dxdydt with an additional 
assumption 

<t:S; (24) 

It is clear that this condition is normal for the collider 
experiments but it is often not so bad for fixed target ex¬ 
periments as well. Due to the smallness of Vm we can keep 
only leading term in the expansion of integrand (0) over 
V. First non-vanishing term contents the first derivative 
of the integrand. The main contribution to this derivative 
comes from e:xjp{byt) (by is a slope parameter). It is not dif¬ 
ficult to keep all the contributions to the derivative but for 
simplicity we restrict our calculation to this contribution. 
In this case 


0^^ 


- I 


' l + T - yi 


(25) 


and it is possible to integrate over (ph, Pk and r explicitly: 


— f yd, ~ + ml 

2ttJ + 

0 

2'7r Tmax 

j dPk J dTF°ji = -2{U - 1). (26) 

0 T-min 

Final result for dp (see eq.(^) is simple 

2oibyVm ,, .. N Q d-my _ ^'o'7^ 

crp = --- (Im - l)o _r, 2 _^ 2 ^ 0 - (27) 

TT Ox ^ ''^v 


Apart from a simple analytical form the obtained for¬ 
mula has one more advantage. The correction depends 
only on the kinematical variables but not on the dynamical 
characteristics of the interaction like ap and ap. It allows 
to eliminate a possible systematical error coming from the 
choice of the model. Such systematics could be large be¬ 
cause we have to know the differential cross sections ap 
and ap in the wide region of varying of four variables, but 
there are neither enough experimental data for that nor 
satisfactory models. Using (|^ ) we have systematics com¬ 
ing from the using of the approximate formulae instead of 
exact ones. But we are able to control it by comparing the 
values calculated by (HI) and (p^ ) using any model. 

Moreover there is a possibility to provide the procedure 
’event by event’ reweighting each event with the RC factor. 
The comparison of the results for the cross sections with 
and without this reweighting gives the correction. 


4 Numerical analysis and code DIFFRAD 

In this section we present the FORTRAN code DIFFRAD 
created on the basis of the exact formulae performed in 


Section ||. The program calculates the lowest order RC to 
the diffractive vector meson electroproduction. The higher 
order effects are approximated by the procedure of expo¬ 
nentiation. The formulae for the cross section are given 
in a covariant form, so the code can be run both for the 
fixed target experiments and for the experiments at col¬ 
lider. The model for ar, t presented originally in ref, 
and developed in ref.||I^ is used as an input. 

Below we give numerical results for RC factor 


_ ^obs 

^ ^0 


277 

J" d4^h^obs 

0 _ 

277 

/ dphC^o 

0 


(28) 


and t- and ^-distributions obtained within the kinematical 
regions of recent experiments on the electro- and muon- 
production of vector mesons. 


4.1 RC to the cross section 

On Figure || one can see Q^-and t-dependences of RC fac¬ 
tor rj in the kinematical region of experiments of EMC 
and NMC. There is no cut on inelasticity used. For high 
values of t the RC can reach factor two. The reason of 
so large effect is smallness of the Born photoproduction 
cross sections aL,T- They fall as exp(byt). In the observed 
cross section this factor is under integral and there is a 
contribution from the region on t where are not so 
small. As a result Gobs falls with the increasing —t but not 
so fast. 

Figure ^ gives the results for rj within kinematics of 
the experiment E665 with a cut on the inelasticity. No¬ 
tice it can be done because of the rather good resolution 
over V (standard deviation a of its distribution is more 
smaller than Vm given by kinematics). The usage of this 
cut leads to a different behavior of ry as a function of Q^. 



g2,GeV2 

Fig. 3. RC factor in kinematical conditions of 
EMC/NMC (pLp^ pLpp). 
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= lOOGeV^ = 200GeV2 


■n 



Q^GeV^ Q^.GeV^ 

Fig. 4. RG factor in kinematical conditions of E665 {fip Fig. 6. RG factor in kinematical conditions of 

ppp). Curves from top to bottom correspond to t =-0.9, -0.5, HERMES {ep epp). 
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Fig. 5. RG factor in kinematical conditions of collider exper¬ 
iments at HERA {ep —> epp). Symbols from top to bottom 
correspond to t =-0.7, -0.5, -0.3, -0.1 GeV^ 


(T, mb 



Fig. 7. The Born (dashed) and observed (solid) cross sections 
in respect to —t. No cut on inelasticity is used. 


The different plots on this figure give IT^-dependence of 

■n- 

The dependence of rj on the kinematical variables 
W'^ and t in the region of the collider experiments at 
HERA is presented in Figure H. No cut on inelasticity 
is used, so Q^-dependence is similar to one given for the 
EMC/NMC experiment. It is found that the RC factor rj 
is not sensitive for . It is due to the fact that the photo¬ 
production cross section is almost flat in the kinematical 
region of the collider experiments. 

The dependence on the Vm cut in the region of HER¬ 
MES kinematics is analyzed in Figure The usage of the 


cut changes RC factor for small and does not influence 
on the RC for larger values of . It can be understood 
from Figure In the case of the cut usage we have to 
define Vm as a minimum of the value of the cut and Vm 
given by the kinematical restrictions (ill)- For fixed values 
of —t the cut influence on w™ and RC till certain value of 
only. 

4.2 t—distribution 

The radiative effects can influence on the slope of the ob¬ 
served cross section in respect to —t due to essential de¬ 
pendence of RC on this variable. It is illustrated in Figures 
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Fig. 8. Slope parameter in respect to —t 

and The t-dependence of the Born cross section within 
the model of refs. 00 is basically defined by two-gluon 
formfactor taken in exponential form exp{byt) with by = 5. 
The dependence on pt also influences on the slope of the 
Born cross section in respect to —t. An additional depen¬ 
dence on t in the case of RC comes basically from Vm, 
that is proportional to f — tmin in the diffractive region of 
small —t. Due to this fact both exponent exp(i5m/) and 
the cross section up tend to zero when t tmin- So the 
observed cross section vanishes in this limit. 


4.3 Inelasticity distribution 

There are several reasons for the inelasticity to be non¬ 
zero: the finite experimental resolution, the mixture of 
non-exclusive events and radiative effects. In this paper 
we assume that a special procedure of separation of exclu¬ 
sive and non-exclusive events has been performed. Strictly 
speaking such procedure can suppress radiative effects par¬ 
tially. It depends on the details of data processing of a 
concrete experiment, and we assume for simplicity that 
there are no such suppression. 

Figure ^ gives two distribution: Born (pure Gaussian 
with a zeroth mean value coming from the effects of a finite 
experimental resolution) and with taking into account RC. 
The second distribution consists of events of two types. 
First ones come from the Born, loop and radiated events 
with V < a. This distribution have the same mean value 
(v = 0) and a as the Born distribution. Other part of 
the events in the radiative corrected distribution comes 
from the radiated events with v > a. Total number of two 
distributions is the same. 


Fig. 9. Inelasticity distributions: Born (line) and radiatively 
corrected (circles) 

5 Discussion and Conclusion 

In this article the QED radiative effects have been an¬ 
alyzed in kinematics of recent experiments on the ex¬ 
clusive vector meson electroproduction. The explicit co¬ 
variant formulae for RC to the cross section are given in 
eqs.d^HI). An approximate expression for the bremsst- 
rahlung cross section can be found in eg. (p^ . 

The RC to the cross section of the diffractive vector 
meson electroproduction is very sensitive to the cut on the 
inelasticity. Using harder cut leads to smaller values of RC 
factor. 

In the diffractive region {—t < 0.3) RC is negative and 
can reach 

— 10% for muon experiments with fixed target; 

— 20% for electron experiments with fixed target; 

— 40% for electron collider experiments. 

This large effect comes basically from a double logarithmic 
contribution in 6inf (0)- For example for collider kine¬ 
matics both its logarithms exceed 10, and 6inf can reach 

0.5. 

There are no essential dependence on the type of a 
vector meson. All numerical results are given for the case 
of p meson production. The dependence on the type of a 
scattered lepton is typical. RC in the case of the electron 
scattering is several times more due to the appearing the 
lepton mass in the argument of the leading logarithm. 

The observed cross section has steeper slope in respect 
to t than the Born cross section. RC to the slope parameter 
is negative and ~ 10%. 

FORTRAN code DIFFRAD is available (aku@hep.by) 
for the calculation of RC to observable quantities in exper¬ 
iments on the diffractive vector meson electroproduction. 
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